Quantum Statistics with Classical Particles 
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Indistinguishability of particles is normally considered to be an inherently quantum property which 
cannot be possessed by a classical theory. However, Saunders Q has argued that this is incorrect, and 
that classically indistinguishable particles are possible. I make this suggestion concrete by describing 
a class of microscopic classical theories involving indistinguishable particles hopping stochastically 
on a graph, and show that it should be possible to experimentally create a physical system realizing 
a simple model by continuously observing atoms trapped in an optical lattice. The indistinguishable 
classical particles obey Bose-Einstein statistics, display the associated clustering phenomena, and 
in appropriate models, can even undergo Bose-Einstein condensation. 



Since the discovery of quantum mechanics, much has 
been made of phenomena such as lasers, Bose-Einstein 
condensation, superfluidity, and superconductivity. The- 
oretical and experimental studies of these phenomena 
have revealed a variety of fascinating and counter- 
intuitive behaviors, consequences of the deviations from 
classical Maxwell-Boltzman statistics displayed by quan- 
tum particles. Indeed, they are associated so strongly 
with quantum mechanics that they are usually described 
as "macroscopic quantum phenomena." Indistinguisha- 
bility, we are told, is a property that can only be pos- 
sessed by quantum particles and not by classical ones. 

Given the importance of these macroscopic quantum 
behaviors, it is worthwhile examining this belief carefully. 
While all our known fundamental quantum particles are 
either bosons or fermions, certainly it is possible in prin- 
ciple to have quantum particles which are distinguish- 
able, perhaps because of some internal quantum number 
which can take a large variety of different values. But is 
it possible to have classical bosons or fermions? 

The standard folklore of physics says "no." Bach 0, 
for instance, has argued that a theory with indistinguish- 
able classical particles is impossible because the trajec- 
tories would distinguish those particles. However, Saun- 
ders has recently argued that one can imagine classical 
indistinguishable particles. Rather than having a sepa- 
rate configuration space for each particle, the space of 
configurations would be a collective one, indicating only 
the number of particles present at each location. In Saun- 
ders' models, these indistinguishable classical particles 
still effectively obeyed Maxwell-Boltzman statistics. The 
difference with quantum particles, he says, is that quan- 
tum particles also experience a concentration of measure 
to a lattice of points (in a system of finite size), whereas 
classical particles can be anywhere. 

However, classical particles which are hopping from 
point to point on a discrete lattice do have their probabil- 
ity measures concentrated on points. In this paper, I give 
explicit microscopic stochastic models with indistinguish- 
able classical particles hopping on a discrete lattice. The 
particles are classical in the sense that at all times, they 
have definite positions. There are no superpositions and 



interference cannot occur. These models display Bose- 
Einstein statistics, unlike Saunders' continuous models, 
and can even display Bose-Einstein condensation, nor- 
mally considered a purely quantum phenomenon. Saun- 
ders also did not specify the microscopic dynamics of his 
particles, leaving open the possibility that a contradiction 
arises somewhere, whereas my models have a complete 
microscopic description. 

The difference between classical indistinguishable par- 
ticles in discrete space and continuous space can be easily 
understood — indistinguishability only has an effect on 
the equilibrium state of a system when particles have a 
good chance to be in the same state. For classical parti- 
cles in a continuous space, the probability is of two par- 
ticles having the same position and momentum, whereas 
quantum particles have a finite probability of overlapping 
due to the uncertainty in position and momentum. On 
a discrete configuration space, however, classical parti- 
cles can arrive at the same state with finite probability, 
and thus can display the distinctive behaviors normally 
associated with quantum indistinguishability. 

Furthermore, the models that I present of classical in- 
distinguishable particles are not just theoretical fantasies; 
it should be possible to realize some simple instances ex- 
perimentally. In particular, spin-0 atoms confined to an 
optical lattice satisfy the Bose-Hubbard model to a good 
approximation 0. This produces a system of bosons 
hopping in a discrete configuration space. These parti- 
cles are quantum, of course, not classical, but suppose 
we observe the system at intervals of length At to see 
how many atoms are present at each lattice site. At 
each observation, the system collapses to some definite 
classical configuration. The required measurement non- 
destructively localizes all particles to individual wells in 
the optical lattice. This, unfortunately, is beyond cur- 
rent experimental techniques, but may not be for long. 
Between observations, some particles may hop from one 
site to another, but if we observe frequently enough, it is 
extremely unlikely that more than one particle will hop 
in time At. Interference between different paths is then 
impossible, and the system becomes, to an excellent ap- 
proximation, a system of classical bosons. 
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After I completed the research reported here, I discov- 
ered similar work by Kaniadakis and Quarati [j, Q • Owe 
models are related but distinct. Their models are primar- 
ily built in continuous systems with an explicit smearing 
function in space and velocity to allow the classical par- 
ticles to overlap (although they do also give one example 
on a lattice). Another similar approach is the quantum 
Boltzmann master equation by Jaksch et al. pj, which 
was used to numerically simulate Bose-Einstein conden- 
sation. I work exclusively in discrete systems, allowing 
my models to be simpler than either of these and to more 
clearly separate quantum mechanical behavior from in- 
distinguishability; indeed, Kaniadakis and Quarati de- 
scribe indistinguishability in their models as a residual 
quantum effect rather than a separate phenomenon, and 
Jaksch et al. make no comment at all about classicality. 

Let us be concrete. Suppose we have some graph G. 
We can associate to it a simple model of classical bosons 
by defining the state of the system to be a function 
n : G — *■ N which gives the number of particles at 
each site i in the graph. The total number of particles 
is TV = Ui. These particles can hop from site to site, 
conserving the total number. For simplicity, we impose 
artificially the constraint that only one particle in the 
whole system can hop in a single time step. This means 
the particles will, strictly speaking, be interacting, but if 
we work in the limit where hopping is a rare event, then 
the interaction is negligible. 

Suppose i and j are two adjacent sites. Then let us 
assume there is a probability p to go in one time step 
from a configuration (1^, 0^ ) (with 1 particle at site i and 

particles at site j) to (Oi, Ij) (with particles at i and 

1 particle at j). This probability is independent of i and 
j (provided they are adjacent) and is the same regardless 
of the number of particles at all other locations. 

However, it is not immediately clear from this what 
the probability should be of hopping from the configura- 
tion {ci,dj) to ((c — {d + l)j). One constraint that 
seems sensible is to require that the probability of hop- 
ping from configuration A to configuration B is the same 
as the probability of hopping from configuration B to 
configuration A. The analogous system of distinguish- 
able classical particles hopping on the graph G has this 
property, since the probability of a particular particle 
hopping forward is the same as its probability of hop- 
ping backwards. This property also corresponds more or 
less to having a Hermitian Hamiltonian in the analogous 
quantum system. 

If we add the condition that the probability of hopping 
out of site i is proportional to the number of particles at 
i, we get the probability rule 

P [ic,,dj) ^ ((c - 1)„ {d + 1),)] = c{d + l)p, (1) 

again assuming the lattice sites i and j are adjacent. We 
should add the additional constraint that p is small, so 
that the total probability of hopping is always less than 



one. Then for any configuration A, we can calculate 
the probability of staying at the same configuration as 
I — J2b -^[^ where the sum is taken over all con- 

figurations B which are one hop away from A. With N 
particles, it is sufficient to take [N^ /4:+{g + l/2)N]p < I, 
where g is the maximum degree of any node in G. Sys- 
tems with distinguishable particles can also obey this 
transition rule, leading to similar effects (e.g., 0)- 

This prescription has the virtue that it agrees with the 
system produced when we frequently observe particles 
obeying a Bose-Hubbard model, as with spin-0 atoms in 
an optical lattice. The standard optical lattice is a cubic 
lattice, but it is also possible to make other graphs by 
blocking some sites using a Fermi-Bose mixture jSJ. The 
Bose-Hubbard model has the Hamiltonian 

H = ^jj2 + J2 + J2 ^''('■~^' - 1)' (2) 

where hi is the annihilation operator at site i and fii — 
h\hi is the number operator at site i. The basis states of 
this quantum system are the possible configurations of 
the above model of classical bosons with the same graph 
of possible sites. If we observe the system at two times 
separated by an interval of At, the probability of hopping 
from a basis configuration A to i? is 



P{A^ B) = {A\e 



JH At/h 



\B) 



(3) 



When At is small, we can approximate 

e'"'^'/" = l + iHAt/h+0{{At)'^). (4) 

The fii terms do not cause transitions, so for A ^ B, 
we only need to consider the term — j) ^l^j 
If A and B differ on the adjacent sites i and j, with A 
including (a, dj) and B including ((c— l)i, (d+ l)j), then 



P{A ^ B) = 



J At 

h 

+0{{Atf 



[c,,dj\b\bj\{c-lhAd+l),) 



(5) 



Recalling that hi\ci) — Vc|(c — 1),), we find 

P{A ^B)= c{d + l){JAt/hf + 0{{Atf), (6) 

which agrees with eq. ^ using p = {J At/h)'^. This ap- 
proximation is valid when gN J At/h, {N At/h) max(ei), 
and N'^U At/h are all much less than 1. 

The classical boson model is a Markov processes, and 
can thus be described via a transition matrix M giving 
the probability to hop between configurations. Because 
the probability of hopping forward is the same as the 
probability of hopping backwards, the transition matrix 
is symmetric, and is thus doubly stochastic (each row 
and each column of M sums to I). For an arbitrary 
initial state, a Markov process asymptotes to an equi- 
librium state, an eigenvector of M with eigenvalue -l-I. 
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For a doubly stochastic matrix, the uniform mixture over 
all configurations is a possible equilibrium state. If the 
graph G is connected, the Markov chain is irreducible, 
and lip is small, there is a non-zero probability of staying 
in the same configuration, so the Markov chain is aperi- 
odic. When both of these are true, the uniform mixture 
will be the only equilibrium state [lOj . 

The equilibrium state of the Markov process there- 
fore corresponds to thermodynamic equilibrium. For in- 
stance, suppose we consider an extremely simple case 
with two lattice sites and two particles. For classical 
distinguishable particles, there are four configurations, 
two with one particle on each site, and two with both 
particles on the same site. The equilibrium state thus 
has probability 1 /2 of having both particles in the same 
location. For the classical boson model, there are only 
three configurations, the same two with both particles 
on the same site, but only one with one particle on each 
site. The equilibrium state thus has probability 2/3 of 
having both particles in the same location, displaying 
the usual clustering effect associated with Bose-Einstein 
statistics. Naturally, this clustering becomes more preva- 
lent in larger systems with a high density of particles. 

This simple model of classical bosons illustrates the ba- 
sic point — that classical particles can be indistinguish- 
able — but cannot display the more exciting properties 
associated with indistinguishable quantum particles. If 
we want to make classical analogues of macroscopic quan- 
tum phenomena, we will need models where the particles 
have additional properties. As a first step, let us intro- 
duce energy into the model. Since we wish energy to 
be conserved at each transition, the new model will have 
two types of particles, which I will call "atoms" and "pho- 
tons." This enables us to have both interesting dynamics 
and conservation of energy by shifting energy back and 
forth between the atoms and the photons. 

Atoms will hop on a graph, just as in the simple model. 
However, now each site i in the graph will have an integral 
value of energy Ei associated with it. We assume that 
adjacent sites of the graph have energy values that are 
within ±1 of each other. (The energies of adjacent sites 
may be the same.) The number of atoms is conserved in 
each transition. We will assume the atoms are classical 
bosons, so the possible atom configurations are described 
by the number of atoms at each site. 

Photons have a graph of their own, and may hop on it 
also as in the previous model. The photons will be indis- 
tinguishable classical particles too. When photons hop 
between different sites on the photon graph, the number 
of photons is conserved, but photons can be created or 
destroyed when an atom makes a transition. In particu- 
lar, we associate a photon site k with each edge (i, j) of 
the atom graph for which \Ei — Ej\ = \. When an atom 
makes a transition along the edge a photon must 

be created or destroyed at the site k in order to conserve 
energy. Each photon is considered to have unit energy. 



Transition type 



Probability 



i ^ j, Ei = Ej niijij + l)p 

i j, Ei — Ej + 1, create 7 at A; ni(nj -f l){mk + l)p' 
i j, Ei — Ej — 1, destroy 7 at fc ni{nj + l)mkp' 



■j, \E,-E,\ > 1 





nikimi + l)q 



TABLE L Transition probabilities in the classical boson en- 
ergy model. Sites i and j are adjacent atom sites with energies 
Ei and Ej. Sites k and I are adjacent photon (7) sites, and 
if appropriate, photon site k is associated to the edge 
There are Ui and rij atoms at sites i and j respectively, and 
nik and mi photons at sites k and I. 



We thus have a variety of possible types of transitions 
in this model, with transition probabilities given in ta- 
ble As with the previous model, we only allow a single 
particle to make a hopping transition in each time step, 
plus the possibility that a single photon is also created 
or destroyed. To determine the transition probabilities, 
it seems reasonable to assume the probability of absorb- 
ing a photon increases proportionally to the number of 
photons available to be absorbed. From this, it follows 
that the requirement that the Markov chain be doubly 
stochastic implies that the probability to emit a pho- 
ton into a site with photons must be proportional to 
TOfc -I- 1. We allow different probabilities p and p' for an 
atom hopping to a site with the same or different energy, 
and a separate probability q for photons hopping. 

Assuming the atom and photon graphs are connected, 
and that p, p', and q are sufficiently small, there is again a 
unique equilibrium state, which is a uniform mixture over 
all configurations. There are two conserved quantities, 
the total energy E and the total number of atoms Na- 
If there are Ui atoms at site i, then the total number 
of photons must be np = E — '^^niEi. Clearly, the 
allowed configurations must satisfy 'Y^^riiEi < E as well 
as ^rti = Nji. Assume there are a total of V photon 
sites. 

We can thus write down a formula for (n^), the ex- 
pected number of atoms at site i. There are 
total possible configurations for np photons, so the total 
number of configurations of all kinds is 



T 



E 



np + V 
np 



1 



(7) 



where the sum is taken over atom configurations {ui} 
satisfying the constraints on total energy and total atom 
number. We thus have 



np + V 
np 



1 



(8) 



To proceed further, we wish to treat the photon system 
as a thermal bath, so we work in the limit where E and V 
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are both very large, but EjV is constant. Then we expect 
most of the energy of the system to be in the photons, so 
the constraint Ea = '^i^i ^ ^ becomes unimportant, 
and we can neglect terms of order Ea/{E + V). Letting 
X — np/{np + V — 1) and xq — E/{E + V), we have 

^'^("""tl ~ i^P + ^Mx) (9) 
w {E + V)h{xo) + EAlnxo,{lO) 

where we have used Stirling's formula to approximate the 
binomial coefficient, andh{x) — — a;lna;— (1 — x) ln(l — x). 
By identifying /3 = — In ccq , we thus produce the standard 
canonical ensemble for the atoms. Indeed, the proce- 
dure above simply mimics the usual argument moving 
from the microcanonical ensemble to the canonical en- 
semble. Since the sum over configurations only considers 
the number of atoms at each site, the atoms behave as 
bosons, not as Maxwell-Boltzman particles. 

We can thus immediately apply the standard results 
about thermodynamics of bosons. For instance, appro- 
priate systems will display Bose-Einstein condensation; 
we need only find an assignment of energies to a graph 
which produces the correct number of states of a given 
energy. For instance, suppose we take the positive octant 
of a square lattice in 3 dimensions, and let the energy of 
site (x, y, z) be x + y + z {x,y,z > 0). We then repli- 
cate the density of states for a 3-dimensional harmonic 
trap with appropriate trap frequency, and in this svstem, 
bosons can undergo Bose-Einstein condensation 

The exact nature of the thermal bath is presumably 
not important, and it could consist of regular quantum 
photons or a different species of atoms. While the model 
discussed above would be difficult to replicate experi- 
mentally, given the variety of possibilities, it might be 
possible to find some system which could experimentally 
realize classical Bose-Einstein condensation. 

Many further extensions are possible. One could con- 
sider interacting particles and could add more parti- 
cle types and more properties such as charge to try to 
replicate other macroscopic quantum phenomena such 
as superconductivity. One could add velocities to make 
bosonic versions of lattice gas automata However, 
deterministic classical boson models face a special chal- 
lenge. In a deterministic model with distinguishable par- 
ticles, two particles with the same properties that start 
with the same state must retain the same state through- 
out their evolution; if they were to separate, there would 
be no way to determine which of the two particles would 
go in which direction. We might worry that this would 
preclude any interesting effects of bosonic particle statis- 
tics, but luckily models with indistinguishable particles 
do not face this problem: It is perfectly possible for two 
indistinguishable particles with the same state to head off 
in different directions without violating any symmetry. 

One could also try to create models of classical 



fermions. A straightforward way of doing this is to sim- 
ply impose a constraint that no site can contain more 
than one of the particles. This prescription immediately 
reproduces the statistical mechanics of fermions using 
the techniques described above for bosons, and as be- 
fore, could be produced experimentally by frequently ob- 
serving quantum fermions moving in an optical lattice. 
However, it is somewhat unsatisfying from a philosophi- 
cal perspective. The constraint causes the models to be- 
have like fermions even if the particles in the model are 
actually distinguishable. It is unclear if there is a true 
distinction between distinguishable and indistinguishable 
classical fermions, but one place to look might be in mod- 
els where two or more particles can hop at the same time 
(as in .41)- Then, for instance, if two particles switch 
locations, this produces a different configuration when 
the particles are distinguishable, but not if they are in- 
distinguishable. Unfortunately, I do not see how such 
models could be realized experimentally, as a double hop 
in a single time step opens the possibility for interference 
between observations of the quantum systems. 

In summary, I have presented some models of classical 
bosons, demonstrating an error in the standard folklore 
that indistinguishability is an inherently quantum prop- 
erty. Indeed, we see that indistinguishability and quan- 
tum behavior are separate phenomena; each can exist 
without the other. In practice, small objects like atoms 
tend to be both quantum and indistinguishable. Larger 
objects have more accessible internal degrees of freedom, 
so tend to lose indistinguishability. Large objects also 
tend to lose quantum coherence, which is perhaps why in- 
distinguishability and quantum behavior have been con- 
sidered in the past to be so closely associated. 

In this paper I have only examined the equilibrium 
behavior of the classical models to show that they can 
reproduce the equilibrium behavior of quantum indistin- 
guishable systems. Of course, it is also possible to study 
the non-equilibrium behavior of classical boson or clas- 
sical fermion models, for instance to examine transport 
properties of the models. Indeed, it is in the realm of dy- 
namics that we can expect to see a difference between the 
quantum and classical models, as interference can play a 
role in the quantum systems but not in the classical ones. 

The classical boson models offer a new perspective for 
understanding macroscopic quantum phenomena. They 
may even provide an arena to make improved concrete 
predictions about such phenomena: Classical systems are 
much easier to simulate and analyze than quantum sys- 
tems, so if classical boson (or fermion) models can be 
created which replicate the major properties of interest- 
ing systems such as high- Tc superconductors, they would 
provide a very useful technique for understanding those 
systems. Indeed, since collective phenomena can persist 
for systems large enough to experience significant deco- 
herence, it might even be that classical boson or fermion 
models are more accurate than existing quantum models. 
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